Quality control and evaluation of plant epigenomics data.

Epigenomics is the study of molecular signatures associated with discrete regions within genomes, many of which are important for a wide range of nuclear processes. The ability to profile the epigenomic landscape associated with genes, repetitive regions, transposons, transcription, differential expression, cis-regulatory elements, and 3D chromatin interactions has vastly improved our understanding of plant genomes. However, many epigenomic and single-cell genomic assays are challenging to perform in plants, leading to a wide range of data quality issues; thus, the data require rigorous evaluation prior to downstream analyses and interpretation. In this commentary, we provide considerations for the evaluation of plant epigenomics and single-cell genomics data quality with the aim of improving the quality and utility of studies using those data across diverse plant species.

RNA-Seq analysis of genes affected by Cyclophilin A/DIAGEOTROPICA (DGT) in tomato root development.

Cyclophilin A/DIAGEOTROPICA (DGT) has been linked to auxin-regulated development in tomato and appears to affect multiple developmental pathways. Loss of DGT function results in a pleiotropic phenotype that is strongest in the roots, including shortened roots with no lateral branching. Here, we present an RNA-Seq dataset comparing the gene expression profiles of wildtype ('Ailsa Craig') and dgt tissues from three spatially separated developmental stages of the tomato root tip, with three replicates for each tissue and genotype. We also identify differentially expressed genes, provide an initial comparison of genes affected in each genotype and tissue, and provide the pipeline used to analyze the data. Further analysis of this dataset can be used to gain insight into the effects of DGT on various root developmental pathways in tomato.

Metabolomics analysis reveals both plant variety and choice of hormone treatment modulate vinca alkaloid production in Catharanthus roseus.

The medicinal plant Catharanthus roseus produces numerous secondary metabolites of interest for the treatment of many diseases - most notably for the terpene indole alkaloid (TIA) vinblastine, which is used in the treatment of leukemia and Hodgkin's lymphoma. Historically, methyl jasmonate (MeJA) has been used to induce TIA production, but in the past, this has only been investigated in whole seedlings, cell culture, or hairy root culture. This study examines the effects of the phytohormones MeJA and ethylene on the induction of TIA biosynthesis and accumulation in the shoots and roots of 8-day-old seedlings of two varieties of C. roseus. Using LCMS and RT-qPCR, we demonstrate the importance of variety selection, as we observe markedly different induction patterns of important TIA precursor compounds. Additionally, both phytohormone choice and concentration have significant effects on TIA biosynthesis. Finally, our study suggests that several early-induction pathway steps as well as pathway-specific genes are likely to be transcriptionally regulated. Our findings highlight the need for a complete set of'omics resources in commonly used C. roseus varieties and the need for caution when extrapolating results from one cultivar to another.

Identification of transcription factors from NF-Y, NAC, and SPL families responding to osmotic stress in multiple tomato varieties.

Identifying osmotic stress-responsive transcription factors (TFs) can facilitate discovery of master regulators mediating salt and/or drought tolerance. To date, few RNA-seq datasets for high resolution time course of salt or drought stress treatments are publicly available for certain crop species. However, such datasets may be available for other crops, and in combination with orthology analysis may be used to infer candidate osmotic stress regulators across distantly related species. Here, we demonstrate the utility of this approach for identification and validation of osmotic stress-responsive transcription factors in tomato. First, we developed physiologically calibrated salt and dehydration-responsive systems for tomato cultivars using real time measurements of transpiration rate and photosynthetic efficiency. Next, we identified differentially expressed TFs in rice using raw RNA-seq datasets for a publicly available salt stress time course. Putative salt stress-responsive TFs in tomato were then inferred based on their orthology with the transcription factors upregulated by salt in rice. Finally, using our osmotic stress system, we experimentally validated stress-responsive expression of predicted tomato candidates representing NUCLEAR FACTOR Y, SQUAMOSA PROMOTER BINDING, and NAC domain TF families. Quantification of transcript copy numbers confirmed that mRNAs encoding all three TFs were strongly upregulated not only by salt but also by drought stress. Induction by both salt and dehydration occurred in a temporal manner across diverse tomato cultivars, suggesting that the identified TFs may play important roles in regulating osmotic stress responses.

NanoCAGE-XL: An Approach to High-Confidence Transcription Start Site Sequencing.

Identifying the transcription start sites (TSS) of genes is essential for characterizing promoter regions. Several protocols have been developed to capture the 5' end of transcripts via Cap-Analysis of Gene Expression (CAGE) or linker-ligation strategies such as Paired-End Analysis of Transcription Start Sites (PEAT), but often require large amounts of tissue. More recently, nanoCAGE was developed for sequencing on the Illumina GAIIx to overcome this limitation. In this chapter, we present the nanoCAGE-XL protocol, the first publicly available adaptation of nanoCAGE for sequencing on recent ultra-high-throughput platforms such as Illumina HiSeq-2000. NanoCAGE-XL provides a method for precise transcription start site identification in large eukaryotic genomes, even in cases where input total RNA quantity is very limited.

IndeCut evaluates performance of network motif discovery algorithms.

Motivation: Genomic networks represent a complex map of molecular interactions which are descriptive of the biological processes occurring in living cells. Identifying the small over-represented circuitry patterns in these networks helps generate hypotheses about the functional basis of such complex processes. Network motif discovery is a systematic way of achieving this goal. However, a reliable network motif discovery outcome requires generating random background networks which are the result of a uniform and independent graph sampling method. To date, there has been no method to numerically evaluate whether any network motif discovery algorithm performs as intended on realistically sized datasets-thus it was not possible to assess the validity of resulting network motifs. Results: In this work, we present IndeCut, the first method to date that characterizes network motif finding algorithm performance in terms of uniform sampling on realistically sized networks. We demonstrate that it is critical to use IndeCut prior to running any network motif finder for two reasons. First, IndeCut indicates the number of samples needed for a tool to produce an outcome that is both reproducible and accurate. Second, IndeCut allows users to choose the tool that generates samples in the most independent fashion for their network of interest among many available options. Availability and implementation: The open source software package is available at https://github.com/megrawlab/IndeCut. Contact: megrawm@science.oregonstate.edu or david.koslicki@math.oregonstate.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

DNase I SIM: A Simplified In-Nucleus Method for DNase I Hypersensitive Site Sequencing.

Identifying cis-regulatory elements is critical in understanding the direct and indirect interactions that occur within gene regulatory networks. Current approaches include DNase-seq, a technique that combines sensitivity to the nonspecific endonuclease DNase I with high-throughput sequencing to identify regions of regulatory DNA on a genome-wide scale. Yet, challenges still remain in processing recalcitrant tissues that have low DNA content. Here, we describe DNase I SIM (for Simplified In-nucleus Method), a protocol that simplifies and facilitates generation of DNase-seq libraries from plant tissues for high-resolution mapping of DNase I hypersensitive sites. By removing steps requiring the use of gel agarose plugs in DNase-seq, DNase I SIM reduces the time required to perform the protocol by at least 2 days, while also making possible the processing of difficult plant tissues including plant roots.

Establishment of Expression in the SHORTROOT-SCARECROW Transcriptional Cascade through Opposing Activities of Both Activators and Repressors.

Tissue-specific gene expression is often thought to arise from spatially restricted transcriptional cascades. However, it is unclear how expression is established at the top of these cascades in the absence of pre-existing specificity. We generated a transcriptional network to explore how transcription factor expression is established in the Arabidopsis thaliana root ground tissue. Regulators of the SHORTROOT-SCARECROW transcriptional cascade were validated in planta. At the top of this cascade, we identified both activators and repressors of SHORTROOT. The aggregate spatial expression of these regulators is not sufficient to predict transcriptional specificity. Instead, modeling, transcriptional reporters, and synthetic promoters support a mechanism whereby expression at the top of the SHORTROOT-SCARECROW cascade is established through opposing activities of activators and repressors.

Small Genetic Circuits and MicroRNAs: Big Players in Polymerase II Transcriptional Control in Plants.

RNA Polymerase II (Pol II) regulatory cascades involving transcription factors (TFs) and their targets orchestrate the genetic circuitry of every eukaryotic organism. In order to understand how these cascades function, they can be dissected into small genetic networks, each containing just a few Pol II transcribed genes, that generate specific signal-processing outcomes. Small RNA regulatory circuits involve direct regulation of a small RNA by a TF and/or direct regulation of a TF by a small RNA and have been shown to play unique roles in many organisms. Here, we will focus on small RNA regulatory circuits containing Pol II transcribed microRNAs (miRNAs). While the role of miRNA-containing regulatory circuits as modular building blocks for the function of complex networks has long been on the forefront of studies in the animal kingdom, plant studies are poised to take a lead role in this area because of their advantages in probing transcriptional and posttranscriptional control of Pol II genes. The relative simplicity of tissue- and cell-type organization, miRNA targeting, and genomic structure make the Arabidopsis thaliana plant model uniquely amenable for small RNA regulatory circuit studies in a multicellular organism. In this Review, we cover analysis, tools, and validation methods for probing the component interactions in miRNA-containing regulatory circuits. We then review the important roles that plant miRNAs are playing in these circuits and summarize methods for the identification of small genetic circuits that strongly influence plant function. We conclude by noting areas of opportunity where new plant studies are imminently needed.

Improved DNase-seq protocol facilitates high resolution mapping of DNase I hypersensitive sites in roots in Arabidopsis thaliana.

BACKGROUND: Identifying cis-regulatory elements is critical in understanding the direct and indirect regulatory mechanisms of gene expression. Current approaches include DNase-seq, a technique that combines sensitivity to the nonspecific endonuclease DNase I with high throughput sequencing to identify regions of regulatory DNA on a genome-wide scale. While this method was originally developed for human cell lines, later adaptations made the processing of plant tissues possible. Challenges still remain in processing recalcitrant tissues that have low DNA content. RESULTS: By removing steps requiring the use of gel agarose plugs in DNase-seq, we were able to significantly reduce the time required to perform the protocol by at least 2 days, while also making possible the processing of difficult plant tissues. We refer to this simplified protocol as DNase I SIM (for simplified in-nucleus method). We were able to successfully create DNase-seq libraries for both leaf and root tissues in Arabidopsis using DNase I SIM. CONCLUSION: This protocol simplifies and facilitates generation of DNase-seq libraries from plant tissues for high resolution mapping of DNase I hypersensitive sites.

NanoCAGE-XL and CapFilter: an approach to genome wide identification of high confidence transcription start sites.

BACKGROUND: Identifying the transcription start sites (TSS) of genes is essential for characterizing promoter regions. Several protocols have been developed to capture the 5' end of transcripts via Cap Analysis of Gene Expression (CAGE) or linker-ligation strategies such as Paired-End Analysis of Transcription Start Sites (PEAT), but often require large amounts of tissue. More recently, nanoCAGE was developed for sequencing on the Illumina GAIIx to overcome these difficulties. RESULTS: Here we present the first publicly available adaptation of nanoCAGE for sequencing on recent ultra-high throughput platforms such as Illumina HiSeq-2000, and CapFilter, a computational pipeline that greatly increases confidence in TSS identification. We report excellent gene coverage, reproducibility, and precision in transcription start site discovery for samples from Arabidopsis thaliana roots. CONCLUSION: nanoCAGE-XL together with CapFilter allows for genome wide identification of high confidence transcription start sites in large eukaryotic genomes.

TIPR: transcription initiation pattern recognition on a genome scale.

MOTIVATION: The computational identification of gene transcription start sites (TSSs) can provide insights into the regulation and function of genes without performing expensive experiments, particularly in organisms with incomplete annotations. High-resolution general-purpose TSS prediction remains a challenging problem, with little recent progress on the identification and differentiation of TSSs which are arranged in different spatial patterns along the chromosome. RESULTS: In this work, we present the Transcription Initiation Pattern Recognizer (TIPR), a sequence-based machine learning model that identifies TSSs with high accuracy and resolution for multiple spatial distribution patterns along the genome, including broadly distributed TSS patterns that have previously been difficult to characterize. TIPR predicts not only the locations of TSSs but also the expected spatial initiation pattern each TSS will form along the chromosome-a novel capability for TSS prediction algorithms. As spatial initiation patterns are associated with spatiotemporal expression patterns and gene function, this capability has the potential to improve gene annotations and our understanding of the regulation of transcription initiation. The high nucleotide resolution of this model locates TSSs within 10 nucleotides or less on average. AVAILABILITY AND IMPLEMENTATION: Model source code is made available online at http://megraw.cgrb.oregonstate.edu/software/TIPR/. CONTACT: megrawm@science.oregonstate.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Alternative splicing in plants: directing traffic at the crossroads of adaptation and environmental stress.

In recent years, high-throughput sequencing-based analysis of plant transcriptomes has suggested that up to ∼60% of plant gene loci encode alternatively spliced mature transcripts. These studies have also revealed that alternative splicing in plants can be regulated by cell type, developmental stage, the environment, and the circadian clock. Alternative splicing is coupled to RNA surveillance and processing mechanisms, including nonsense mediated decay. Recently, non-protein-coding transcripts have also been shown to undergo alternative splicing. These discoveries collectively describe a robust system of post-transcriptional regulatory feedback loops which influence RNA abundance. In this review, we summarize recent studies describing the specific roles alternative splicing and RNA surveillance play in plant adaptation to environmental stresses and the regulation of the circadian clock.

Environmental stresses modulate abundance and timing of alternatively spliced circadian transcripts in Arabidopsis.

Environmental stresses profoundly altered accumulation of nonsense mRNAs including intron-retaining (IR) transcripts in Arabidopsis. Temporal patterns of stress-induced IR mRNAs were dissected using both oscillating and non-oscillating transcripts. Broad-range thermal cycles triggered a sharp increase in the long IR CCA1 isoforms and altered their phasing to different times of day. Both abiotic and biotic stresses such as drought or Pseudomonas syringae infection induced a similar increase. Thermal stress induced a time delay in accumulation of CCA1 I4Rb transcripts, whereas functional mRNA showed steady oscillations. Our data favor a hypothesis that stress-induced instabilities of the central oscillator can be in part compensated through fluctuations in abundance and out-of-phase oscillations of CCA1 IR transcripts. Taken together, our results support a concept that mRNA abundance can be modulated through altering ratios between functional and nonsense/IR transcripts. SR45 protein specifically bound to the retained CCA1 intron in vitro, suggesting that this splicing factor could be involved in regulation of intron retention. Transcriptomes of nonsense-mediated mRNA decay (NMD)-impaired and heat-stressed plants shared a set of retained introns associated with stress- and defense-inducible transcripts. Constitutive activation of certain stress response networks in an NMD mutant could be linked to disequilibrium between functional and nonsense mRNAs.

The cyclophilin A DIAGEOTROPICA gene affects auxin transport in both root and shoot to control lateral root formation.

Cyclophilin A is a conserved peptidyl-prolyl cis-trans isomerase (PPIase) best known as the cellular receptor of the immunosuppressant cyclosporine A. Despite significant effort, evidence of developmental functions of cyclophilin A in non-plant systems has remained obscure. Mutations in a tomato (Solanum lycopersicum) cyclophilin A ortholog, DIAGEOTROPICA (DGT), have been shown to abolish the organogenesis of lateral roots; however, a mechanistic explanation of the phenotype is lacking. Here, we show that the dgt mutant lacks auxin maxima relevant to priming and specification of lateral root founder cells. DGT is expressed in shoot and root, and localizes to both the nucleus and cytoplasm during lateral root organogenesis. Mutation of ENTIRE/IAA9, a member of the auxin-responsive Aux/IAA protein family of transcriptional repressors, partially restores the inability of dgt to initiate lateral root primordia but not the primordia outgrowth. By comparison, grafting of a wild-type scion restores the process of lateral root formation, consistent with participation of a mobile signal. Antibodies do not detect movement of the DGT protein into the dgt rootstock; however, experiments with radiolabeled auxin and an auxin-specific microelectrode demonstrate abnormal auxin fluxes. Functional studies of DGT in heterologous yeast and tobacco-leaf auxin-transport systems demonstrate that DGT negatively regulates PIN-FORMED (PIN) auxin efflux transporters by affecting their plasma membrane localization. Studies in tomato support complex effects of the dgt mutation on PIN expression level, expression domain and plasma membrane localization. Our data demonstrate that DGT regulates auxin transport in lateral root formation.

Environmental Stresses Modulate Abundance and Timing of Alternatively Spliced Circadian Transcripts in Arabidopsis.

Environmental stresses profoundly altered accumulation of nonsense mRNAs including intron retaining (IR) transcripts in Arabidopsis. Temporal patterns of stress-induced IR mRNAs were dissected using both oscillating and non-oscillating transcripts. Broad range thermal cycles triggered a sharp increase in the long intron retaining CCA1 isoforms and altered their phasing to different times of day. Both abiotic and biotic stresses such as drought or P. syringae infection induced similar increase. Thermal stress induced a time delay in accumulation of CCA1 I4Rb transcripts whereas functional mRNA showed steady oscillations. Our data favor a hypothesis that stress-induced instabilities of the central oscillator can be in part compensated through fluctuations in abundance and out of phase oscillations of CCA1 IR transcripts. Altogether, our results support a concept that mRNA abundance can be modulated through altering ratios between functional and nonsense/IR transcripts. SR45 protein specifically bound to the retained CCA1 intron in vitro, suggesting that this splicing factor could be involved in regulation of intron retention. Transcriptomes of NMD-impaired and heat-stressed plants shared a set of retained introns associated with stress- and defense-inducible transcripts. Constitutive activation of certain stress response networks in an NMD mutant could be linked to disequilibrium between functional and nonsense mRNAs.

A comparative study of ripening among berries of the grape cluster reveals an altered transcriptional programme and enhanced ripening rate in delayed berries.

Transcriptional studies in relation to fruit ripening generally aim to identify the transcriptional states associated with physiological ripening stages and the transcriptional changes between stages within the ripening programme. In non-climacteric fruits such as grape, all ripening-related genes involved in this programme have not been identified, mainly due to the lack of mutants for comparative transcriptomic studies. A feature in grape cluster ripening (Vitis vinifera cv. Pinot noir), where all berries do not initiate the ripening at the same time, was exploited to study their shifted ripening programmes in parallel. Berries that showed marked ripening state differences in a véraison-stage cluster (ripening onset) ultimately reached similar ripeness states toward maturity, indicating the flexibility of the ripening programme. The expression variance between these véraison-stage berry classes, where 11% of the genes were found to be differentially expressed, was reduced significantly toward maturity, resulting in the synchronization of their transcriptional states. Defined quantitative expression changes (transcriptional distances) not only existed between the véraison transitional stages, but also between the véraison to maturity stages, regardless of the berry class. It was observed that lagging berries complete their transcriptional programme in a shorter time through altered gene expressions and ripening-related hormone dynamics, and enhance the rate of physiological ripening progression. Finally, the reduction in expression variance of genes can identify new genes directly associated with ripening and also assess the relevance of gene activity to the phase of the ripening programme.

Paired-end analysis of transcription start sites in Arabidopsis reveals plant-specific promoter signatures.

Understanding plant gene promoter architecture has long been a challenge due to the lack of relevant large-scale data sets and analysis methods. Here, we present a publicly available, large-scale transcription start site (TSS) data set in plants using a high-resolution method for analysis of 5' ends of mRNA transcripts. Our data set is produced using the paired-end analysis of transcription start sites (PEAT) protocol, providing millions of TSS locations from wild-type Columbia-0 Arabidopsis thaliana whole root samples. Using this data set, we grouped TSS reads into "TSS tag clusters" and categorized clusters into three spatial initiation patterns: narrow peak, broad with peak, and weak peak. We then designed a machine learning model that predicts the presence of TSS tag clusters with outstanding sensitivity and specificity for all three initiation patterns. We used this model to analyze the transcription factor binding site content of promoters exhibiting these initiation patterns. In contrast to the canonical notions of TATA-containing and more broad "TATA-less" promoters, the model shows that, in plants, the vast majority of transcription start sites are TATA free and are defined by a large compendium of known DNA sequence binding elements. We present results on the usage of these elements and provide our Plant PEAT Peaks (3PEAT) model that predicts the presence of TSSs directly from sequence.

Sustained-input switches for transcription factors and microRNAs are central building blocks of eukaryotic gene circuits.

WaRSwap is a randomization algorithm that for the first time provides a practical network motif discovery method for large multi-layer networks, for example those that include transcription factors, microRNAs, and non-regulatory protein coding genes. The algorithm is applicable to systems with tens of thousands of genes, while accounting for critical aspects of biological networks, including self-loops, large hubs, and target rearrangements. We validate WaRSwap on a newly inferred regulatory network from Arabidopsis thaliana, and compare outcomes on published Drosophila and human networks. Specifically, sustained input switches are among the few over-represented circuits across this diverse set of eukaryotes.